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ABSTRACT 

Context. The standard ACDM cosmological model implies that all celestial bodies are embedded in a perfectly uniform 
dark energy background, represented by Einstein's cosmological constant, and experience its repulsive antigravity action. 
Aims. Can dark energy have strong dynamical effects on small cosmic scales as well as globally? Continuing our efforts 
to clarify this question, we focus now on the Virgo Cluster and the flow of expansion around it. 

Methods. We interpret the Hubble diagram, from a new database of velocities and distances of galaxies in the cluster 
and its environment, using a nonlinear analytical model which incorporates the antigravity force in terms of Newtonian 
mechanics. The key parameter is the zero-gravity radius, the distance at which gravity and antigravity are in balance. 
Results. 1. The interplay between the gravity of the cluster and the antigravity of the dark energy background determines 
the kinematical structure of the system and controls its evolution. 2. The gravity dominates the quasi-stationary bound 
cluster, while the antigravity controls the Virgocentric flow, bringing order and regularity to the flow, which reaches 
linearity and the global Hubble rate at distances J> 15 Mpc. 3. The cluster and the flow form a system similar to the 
Local Group and its outflow. In the velocity-distance diagram, the cluster-flow structure reproduces the group-flow 
structure with a scaling factor of about 10; the zero-gravity radius for the cluster system is also 10 times larger. 
Conclusions. The phase and dynamical similarity of the systems on the scales of 1-30 Mpc suggests that a two-component 
pattern may be universal for groups and clusters: a quasi-stationary bound central component and an expanding outflow 
around it, due to the nonlinear gravity-antigravity interplay with the dark energy dominating in the flow component. 
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. 1. Introduction 

Early studies of the motion of our Galaxy (reviewed by 
Huchra 1988) led to the discovery of the peculiar velocity 
of the Local Group towards the Virgo cluster. To illustrate, 
take the distance to Virgo to be 17 Mpc and the Hubble 
constant Ho = 72 km/s/Mpc, then the expected cosmo- 
logical velocity of the cluster is 1224 km/s. The observed 
velocity is about 1000 km/s. The difference reflects our pe- 
culiar velocity of about 220 km/s towards Virgo. 

This estimate of the retarded expansion (the so-called 
"Virgo infall") views the universal expansion as existing 
. not only on "truly cosmological" scales ~ 1000 Mpc, but 
also on local scales of ~ 10 Mpc only weakly disturbed. 
Observational and theoretical aspects of the Virgo infall 
were studied e.g. by Silk (1974, 1977), Peebles (1977), 
Hoffman & Salpeter (1982), Sandage (1986), Teerikorpi et 
al. (1992), and Ekholm et al. (1999, 2000). 

1.1. The Virgocentric flow 

The Virgo infall is a particular feature of what is known 
as the Virgocentric flow: hundreds of galaxies (our Galaxy 
among them) are receeding away from the Virgo cluster. 
Recent data (Karachentsev & Nasonova 2010) show that 
the flow velocities range from nearly zero around 5-8 Mpc 
from the cluster center to about 2 000 km/s at the distances 



of about 30 Mpc. For R > 10—15 Mpc, a roughly linear 
velocity-distance relation may be seen in the flow, making 
it resemble the global Hubble expansion. But strong devia- 
tions from the linear relation exist at smaller distances, not 
surprisingly: 1) the linear expansion flow, V — HR, is the 
property of the universe beyond the cosmic "cell of unifor- 
mity" (> 300 Mpc); 2) the Virgo cluster is a strong local 
overdensity; 3) the matter distribution around the cluster 
is highly nonuniform up to the distance of, say, 30 Mpc. 

Besides the natural deviations, we have here an exam- 
ple of a paradox realized by Sandage (1986, 1999). He saw 
a mystery in the fact that the local expansion proceeds 
in a quite regular way. Moreover, the rate of expansion 
is similar, if not identical, to the global Hubble constant 
(Ekholm et al. 2001; Thim et al. 2003; Karachentsev 2005; 
Karachentsev et al. 2002, 2003b; Whiting 2005). Originally, 
the Friedmann theory described the expansion of a smooth, 
uniform self-gravitating medium. But the linear cosmolog- 
ical expansion was discovered by Hubble where it should 
not be: in the lumpy environment at distances < 20 Mpc. 

1.2. The local relevance of dark energy 

Soon after the discovery of dark energy (Riess et al. 1998; 
Perlmutter et al. 1999) on global scales, we suggested 
(Chernin et al. 2000, Chernin 2001, Baryshev et al. 2001, 
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Karachentsev et al. 2003a, Teerikorpi et al. 2005, Chernin 
et al. 2004) that dark energy with its omnipresent and uni- 
form density (represented by Einstein's cosmological con- 
stant) may provide the dynamical background for a regular 
quiescent expansion on local scales, resolving the Hubble- 
Sandage paradox. Our key argument came from the fact 
that the antigravity produced by dark energy is stronger 
than the gravity of the Local Group at distances larger 
than about 1.5 Mpc from the group center [3 

The present paper extends our work from groups to the 
scale of clusters, focusing on the Virgo cluster. 

2. The phase space structure 

Recent observations of the Virgo cluster and its vicinity 
permit a better understanding of the physics behind the 
visible structure and kinematics of the system. 

2.1. The dataset 

The most complete list of data on the Virgo Cluster and 
the Virgocentric flow has been collected by Karachentsev & 
Nasonova (2010). These include distance moduli of galax- 
ies from the Catalogue of the Neighbouring Galaxies (— 
CNG, Karachentsev et al. 2005) and also from the litera- 
ture with the best measurements prefered. Distances from 
the Tip of the Red Giant Branch (TRGB) and the Cepheids 
are used from the CNG together with new TRGB dis- 
tances (Karachentsev et al. 2006, Tully et al. 2006, Mei 
et al. 2007). For galaxy images in two or more photomet- 
ric bands obtained with WFPC2 or ACS cameras at the 
HST, the TRGB method yields distances with an accuracy 
of about 7% (Rizzi et al. 2007). The database includes also 
data on 300 E and SO galaxies from the Surface Brightness 
Fluctuation (SBF) method by Tonry et al. (2000) with a 
typical distance error of 12%. The total sample contains 
the velocities and distances of 1371 galaxies within 30 Mpc 
from the Virgo cluster center. Especially interesting is the 
sample of 761 galaxies selected to avoid the effect of un- 
known tangential (to the line of sight) velocity components. 
The velocity-distance diagram for this sample taken from 
Karachentsev & Nasonova (2010) is given in Fig.l. 

2.2. The zero-velocity radius and the cluster mass 

The zero- velocity radius within the retarded expansion field 
around a point-like mass concentration means the distance 
where the radial velocity relative to the concentration is 
zero. In the ideal case of the mass concentration at rest 
within the expanding Friedmann universe this is the dis- 
tance where the radial peculiar velocity towards the con- 
centration is equal to the Hubble velocity for the same dis- 
tance. Using Tully-Fisher distances in the Hubble diagram, 
Teerikorpi et al. (1992) could for the first time see the lo- 
cation of the zero-velocity radius Ro for the Virgo system, 

1 The increasing observational evidence and theoretical con- 
siderations favoring this view have been discussed by Maccio 
et al. (2005), Sandage (2006), Teerikorpi et al. (2006, 2008), 
Chernin et al. (2006, 2007a,b,c), Byrd et al. (2007), Valtonen 
et al. (2008), Balaguera-Antolinez et al. (2007), Bambi (2007), 
Chernin (2008), Niemi & Valtonen (2009), Guo & Shan (2009). 
For some counter-arguments to this new approach see Hoffman 
et al. (2008) and Martinez- Valquero et al. (2009). 
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Fig. 1. The Hubble diagram with Virgocentric velocities 
and distances for 761 galaxies of the Virgo Cluster and 
the outflow (Karachentsev & Nasonova 2010). The two- 
component pattern is outlined by bold dashed lines. The 
zero-gravity radius Rzg ~ 10 Mpc is located in the zone 
between the components. The broken line is the running 
median used by K & N to find the zero-velocity radius 
Rq = 6 — 7 Mpc. Two lines from the origo indicate the 
Hubble ratios H — 72 km/s/Mpc and H v = 65 km/s/Mpc. 
Two curves show special trajectories, corresponding to the 
parabolic motion (E = 0; the upper one) and to the min- 
imal escape velocity from the cluster potential well (the 
lower one). The flow galaxies with the most accurate dis- 
tances and velocities (big dots) occupy the area between 
these curves. 



so that Ro/Rvirgo « 0.45 or R w 7.4 Mpc. The work by 
Karachentsev & Nasonova (2010) puts the zero- velocity ra- 
dius at R = 5.0 — 7.5 Mpc. For R < R , positive and 
negative velocities appear in practically equal numbers; for 
R > i?0i the velocities are positive with a few exceptions 
likely due to errors in distances. The zero-velocity radius 
gives the upper limit of the size of the gravitationally bound 
cluster, and the diagram shows that the Virgocentric flow 
starts at R > Rq and extends at least up to 30 Mpc. 

The zero- velocity radius has been often used for estimat- 
ing the total mass Mq of a gravitationally bound system. 
According to Lynden-Bell (1981) and Sandage (1986), the 
spherical model with A = leads to the estimator 

Mo = (ir 2 /8G)tu 2 Rl (1) 

With the age of the universe tjj — 13.7 Gyr (Spergel et 
al. 2007), Karachentsev & Nasonova (2010) find the Virgo 
cluster mass Mq = (6.3 ±2.0) x 10 14 M©. This result agrees 
with the virial mass M v ; r = 6 x 10 14 M Q estimated by 
Hoffman & Salpeter (1982) and - 4 x 10 14 M Q of Valtonen 
et al. (1985) and Saarinen & Valtonen (1985). Teerikorpi 
et al. (1992) and Ekholm et al. (1999, 2000) found that the 
real cluster mass might be from 1 to 2 the virial mass, or 
(0.6 - 1.2) x 10 15 M o . Tully & Mohayaee (2004) derived the 
mass 1.2 x 1O 15 M with the "numerical action" method. 
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2.3. The Virgo infall 

Figure 1 shows clearly the Virgo infall: at distances R > 
Ro, the galaxies are located mostly below the Hubble line 
Vh = HR with H = 72 km/s/Mpc (Spergel et al. 2007). 
The deviation of the median velocity of the flow V m from 
the Hubble velocity is well seen in the distance range 7-13 
Mpc. At larger distances the retardation effect gets weaker 
and gradually sinks below the measurement error level. 

For the distance and peculiar velocity of the Milky Way, 
we take figures from Karachentsev & Nasonova (2010): the 
distance to the Virgo cluster 17 Mpc, the recession veloc- 
ity 1004 ± 70 km/s; with the Hubble constant H = 72 
km/s/Mpc the regular expansion velocity is 1224 km/s. 
Then the peculiar velocity directed to the cluster center 
SV = 220 ± 70 km/s (this agrees with the result based on 
TF distances and a different analysis by Theureau et al. 
1997; a complete velocity field on the scale of 10-20 Mpc is 
given by Tully et al. 2008.) The effect is not very strong. It 
is much more instructive, however, that there are big non- 
linear deviations from the linear velocity-distance law near 
the cluster: the infall effect is strongest there. 



2.4. Two-component phase structure 

The Virgo system reveals an obvious two-component struc- 
ture in Fig.l: the cluster and the flow, outlined roughly 
with bold dashed lines. As we noted, the positive and neg- 
ative velocities (from -2000 to +1700 km/s) are seen in the 
cluster area in equal numbers, so the component is rather 
symmetrical relative to the horizontal line V = 0. The bor- 
der zone between the components, in the range 6 < R < 12 
Mpc, is poorly populated, and the velocities are consider- 
ably less scattered here, from —400 to +600 km/s. This 
zone contains the zero- velocity radius. The flow component 
with positive velocities is at distances > 12 Mpc. It is rather 
symmetrical relative to the line V = HR. The velocities are 
scattered within ±1000 km/s around the symmetry line. 

A much smaller velocity dispersion is seen in the more 
accurate data on 75 galaxies with TRGB and Cepheid dis- 
tances (mostly from CNG) located from 10 to 25 Mpc from 
the cluster center. In the diagram of Fig.l, the subsam- 
ple (dots) occupies a strip less than 500 km/s wide. They 
are scattered with the mean dispersion of about 250 km/s 
around the line V = HyR, where Hy ~ 65 km/s/Mpc. 
This line is the median line for the subsample. 
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Fig. 2. The Hubble diagram for 57 galaxies of the Local 
Group and the local flow (Karachetsev et al. 2009). The 
thin line indicates the zero-velocity radius Ro = 0.7 — 0.9 
Mpc. The two-component phase structure is outlined by the 
bold dashed line. Zoomed in about 10 times, the group- 
scale structure reproduces well the cluster-scale structure 
(Fig.l). The zero-gravity radius Rzg is about 10 times 
smaller for the group than the same key quantity of the 
cluster system. The similarity supports a universal cosmic 
two-component grand-design for groups and clusters. 



factor of 10, would roughly reproduce the phase structure 
of the Virgo system. We will see that the phase similarity 
is also supported by the dynamical similarity. 



3. Gravity-antigravity interplay 

Dark energy is a form of cosmic energy that produces 
antigravity and causes the global cosmological accelera- 
tion discovered by Riess et al. (1998) and Perlmutter et 
al. (1999) in observations of SNe type la at horizon-size 
distances of about 1000 Mpc. These and other observa- 
tions, especially the cosmic microwave background (CMB) 
anisotropies (Spergel et al. 2007), indicate that the global 
dark energy density, p v = (0.73 ± 0.03) x 10~ 29 g/cm 3 , 
makes nearly 3/4 of the total energy of the universe. 



2.5. Similarity with the Local system 

The two-component phase structure of the Virgo sys- 
tem looks very similar to that of the Local system (the 
Local Group and the outflow). The diagram of Fig. 2 
(Karachentsev et al. 2009) is based on the most complete 
and accurate data obtained with the TRGB method using 
the HST. The velocities within the Local Group range from 
— 150 to +160 km/s, about 10 times less than the velocity 
spread in the Virgo cluster. The zone between the two com- 
ponents of the Local system is in the range 0.7 — 1.2 Mpc, 
also 10 times less than in Fig.l. The local flow has a dis- 
persion of about 30 km/s around the line V = H^R, where 
Hl = Hy — 60 km/s/Mpc. This line defines the median for 
the local flow. The dispersion is about 10 times less than 
in the Vigocentric flow for the subsample of the accurate 
distances. The local phase diagram of Fig. 2, if zoomed by a 



3.1. Antigravity in Newtonian description 

According to the simplest, most straightforward and quite 
likely view adopted in the ACDM cosmology, dark energy is 
represented by Einstein's cosmological constant A, so that 
pv = (c 2 /8ttG)A. If so, then dark energy is the energy of 
the cosmic vacuum (Glincr 1965) and is described macro- 
scopically as a perfectly uniform fluid with the equation 
of state pv = —pv (here pv is the dark energy pressure; 
c = 1). This interpretation implies that dark energy exists 
everywhere in space with the same density and pressure. 

Local dynamical effects of dark energy have been 
studied by Chernin et al. (2000, 2009), Chernin (2001), 
Baryshev et al. (2001), Karachentsev et al. (2003), Byrd et 
al. (2007), and Teerikorpi et al. (2008). For the Local Group 
and the outflow close to it, we showed that the antigrav- 
ity produced by the dark energy background exceeds the 
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gravity of the group at distances beyond ~ 1.5 Mpc from 
the group barycentcr. Similar results were obtained for the 
nearby M81 and Cen A groups (Chernin et al. 2007a, b). 

Extending now our approach to the cluster scale, we 
consider the Virgo system (the cluster and the flow) em- 
bedded in the uniform dark energy background. We de- 
scribe the system - in the first approximation - in a frame- 
work of a spherical model (Silk 1974, 1977; Peebles 1976) 
where the cluster is considered as a spherical gravitationally 
bound quasi-stationary system. The galaxy distribution in 
the flow is represented by a continuous dust-like (pressure- 
less) matter. The flow is expanding, and it is assumed that 
the concentric shells of matter do not intersect each other 
in their motion (the mass within each shell keeps constant). 

Since the velocities in the cluster and the flow are much 
less than the speed of light, we may analyze the dynamics 
of the cluster-flow system in terms of Newtonian forces. 
A spherical shell with a Lagrangian coordinate r and a 
Eulerian distance R(r, t) from the cluster center is affected 
by the gravitation of the cluster and the flow galaxies within 
the sphere of the radius R(r 7 1). According to the Newtonian 
gravity law, the shell experiences an acceleration 



GM(r) 
R(r,t) 2 



(2) 



in the reference frame related to the cluster center. Here 
M{r) is the total mass within the radius R(r,t). 

The flow is affected also by the repulsion produced 
by the local dark energy. The resulting acceleration in 
Newtonian terms is given by the " Einstein antigravity law" : 



F E (r,t) = +G2 PA ^0 = + ^Gp A R(r,t) 



(3) 



Here — 2p A = pA + 3p A is the effective (General Relativity) 
gravitating density of dark energy (for details see, e.g., 
Chernin 2001). It is negative, and therefore the accelera- 
tion is positive, speeding up the shell away from the center. 

3.2. Zero-gravity radius in the point mass approximation 

The equation of motion of an individual shell has the form 



R(r,t)=F N +F E 



GM(r) 8tt „, . 
■ W ^ + Y G PA R(r,t). (4) 



Its integration leads to the mechanical energy conservation 
law for each matter shell: 



l v(rt) 2 GM(r) , 4tt^ 
2 V[r,t) ~ R(r,t) ' 



-Gp A R(r,tY + E(r). 



(5) 



Here V(r, t) = R(r, t) is the Eulerian radial velocity of 
the shell, and E(r) is its total mechanical energy (per unit 
mass). This is the basic equation which describes the out- 
flow in our model. It is the Newtonian analogue of the LTB 
spacetime of General Relativity (e.g., Gromov et al. 2001). 

A simpler model results, if we can assume that the clus- 
ter mass Mo makes the majority of the mass M(r), and 
the flow galaxies can be treated as test particles. Curiously 
enough, this assumption makes the model more general in 
the sense that no spherical symmetry is required now for 
the galaxy distribution in the flow (not necessarily concen- 
tric shells around the cluster). However, the gravity force 



produced by the cluster at the distances of the flow is still 
isotropic. The model remains also nonlinear: no restriction 
on the strength of the dynamical effects is assumed. 

The mass-point model may seem overly simple, as it 
is known that the Virgo cluster contributes only 20 % of 
the bright galaxies within the volume up to the distance of 
the Local Group (Tully 1982). However, Tully & Mohayaee 
(2004) conclude from their numerical action method that 
the mass-to-luminosity ratio of the E-galaxy rich Virgo clus- 
ter must be several times higher (they give M/L ~ 900) 
than in the surrounding region (~ 125) in order to pro- 
duce the infall velocities observed close to Virgo (a quali- 
tatively similar result was found by Teerikorpi et al. 1992). 
Therefore, the model with most of the mass in and close to 
the Virgo cluster could be a sufficient first approximation. 

In this simplified/generalized version, the mass M(r) = 
M is constant in Eq. 4, which shows that antigravity and 
gravity balance each other at the distance R = i?zGj where 

R ™ = & 1/3 ^ 
HirpA 

is the zero-gravity radius (Chernin et al. 2000; Chernin 
2001; Dolgachev et al. 2003,2004). Gravity dominates at 
R < Rzg, an d antigravity is stronger than gravity at 
R > Rzg - The zero-gravity radius is the major quantitative 
factor in the gravity-antigravity interplay in the system. 

According to various estimates (Sect.2), the total mass 
of the Virgo cluster is 1 - 2 times the virial mass: M = 
(0.6 - 1.2) x 1O 15 M . Then the known global p A gives 



i? ZG = (9-ll) Mpc. 



(7) 



Karachentsev & Nasonova (2010) argue in favour of one 
virial mass for the cluster mass; if so the lower limit for 
Rzg in Eq.(7) may be more realistic. 

Eqs.6,7 indicate that the zero-gravity radius Rzg of 
the Virgo system is not far from the zero-velocity radius 
i?o = 5.0 — 7.5 Mpc. In fact, one expects that Rzg is 
slightly larger than i? (Teerikorpi et al. 2008) in the stan- 
dard cosmology. The zero- velocity radius constrains the size 
of the cluster (Sec. 2), and there also gravity dominates 
(R < Rq SS Rzg), the necessary condition for a gravitation- 
ally bound cluster. On the other hand, the Virgocentric flow 
is located beyond the distance Rzg — Ro, hence the flow is 
dominated dynamically by the dark energy antigravity. 

The zero-gravity radius of the Virgo system is about 10 
times larger than that of the Local system for which Rzg = 
1-1.3 Mpc (Chernin 2001; Chernin et al. 2000, 2009). We 
see here the same factor 10 as for the phase structures of the 
Virgo and Local systems (Sec. 2). Thus the phase similarity 
is complemented by the dynamical similarity. 

3.3. Antigravity push 

Now we return to the Milky Way and estimate the net 
gravity-antigravity effect at its distance i?MW = 17 Mpc 
from the cluster center. With Rzg = 9 — 11 Mpc, we find 
that here the antigravity force is stronger than the gravity: 



F E /|F N | -4-7, R = R 



MW- 



(8) 



As we see, the extra gravity produced by the Virgo cluster 
mass is actually overbalanced by the antigravity of the dark 
energy located within the sphere of the radius Rmw around 
the cluster. Thus the Milky Way experiences not a weak 
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extra pull (cf. Sees. 1,2), but rather a strong push (in the 
Virgocentric reference frame) from the Virgo direction. 

The dark energy domination in the flow changes the 
current dynamical situation drastically in comparison with 
the earlier models by Silk (1974, 1977) and Peebles (1976). 
At the periphery of the Virgocentric flow (at about 30 Mpc) 
the force ratio is impressively large: Fe/\Fn\ — 20 — 37. 

3.4. Antigravity domination 

The basic trend of the flow evolution is seen from Eq.5. 
With M{r) = M = const, it takes at large times the form: 



V{r,t = ^) = [^Gp A ] 1 ' 2 R. 



(9) 



In this limit, -Fe/I-FnI 00 , an d the gravity of the cluster 
is practically negligible compared to the dark energy anti- 
gravity. As a result, the flow acquires asymptotically the 
linear velocity-distance relation with the constant Hubble 
factor expressed via the dark energy density only: 



r 8nG 



-PA 



|l/2 



63 - 65 km/s/Mpc, 



(10) 



At a given moment of time (that of observation) t = to Eq.5 
presents the spatial structure of the flow. At large distances 



V(r = oo, to) -> H A R(r,t ). 



The flow trajectories converge to the asymptotic of Eq.ll 
independently, in general, of the "boundary /initial condi- 
tions" near the cluster (see figure 2 in Teerikorpi & Chernin 
2010). In particular, strong perturbations at small distances 
in the flow are compatible with its regularity at larger dis- 
tances. No tuning of the model parameters are required for 
this - the only significant quantity is the zero-gravity radius 
i?zc determined by the cluster mass and the DE density. 

Thus, dark energy dominates the dynamics of the 
Virgocentric flow on its entire observed extension from 12 
to at least 30 Mpc. Because of the strong antigravity, the 
flow proves to be rather regular acquiring the linear Hubble 
relation at distances J; 15 Mpc from the cluster center. This 
major conclusion from the model also resolves the Hubblc- 
Sandage paradox on the scale of ~ 10 Mpc, in harmony 
with our general approach to the old puzzle (Sec.l). 

4. Discussion 

We discuss here some specific features and implications of 
the dark energy domination in the Virgocentric flow. 

4.1. The Einstein-Straus radius 

The mass of dark matter and baryons in the Virgo Cluster, 
M = (0.6 - 1.2) x 10 15 M Q , was initially collected from the 
volume which at the present epoch has the radius Res 



3M 



11/3 



U7rp M (to) J — 15 — 19 Mpc, where pm(^o) is the mean 
matter density in the Universe. The quantity Res is known 
as the Einstein-Straus radius in the vacuole model. In the 
presence of dark energy (Chernin et al. 2005; Teerikorpi 
et al. 2008), R ES = [ 1 ^} 1/3 Rzg = 1.7 Rzg- Within the 
E-S vacuole, the spherical layer between the radii Ro and 
Res does not contain any mass. Therefore in this layer, 
the approximation of test particles adopted in the model 



of Sec. 3 works well. Beyond the distance -Res — Rmw, the 
model needs modification, but remains qualitatively correct 
up to the distances of 25-30 Mpc. Anyway since 2p A > 
Pnito), the antigravity of dark energy dominates in the 
entire area of the Virgocentric flow in Fig.l. 

4.2. Two special trajectories 

Eq.5 shows that the flow has zero velocity at R = Rzg — 
Ro, if the total energy of the shell or an individual galaxy 



E(r) = E min = - 



3 GM 



2 Rzc 



(12) 



It is easy to see that E m i n is the minimal energy needed for 
a particle to escape from the potential well of the cluster. 
Eq.5 with E(r) = E m i D gives the corresponding velocity 
V m i n (R) of the flow at different distances in the flow: 



V mi „ = H A R[l + 2( 



R ' 



3(^f)T /2 - 



(13) 



At large distances (R » Rzg), the trajectories with 
E(r) = E m i n approach the asymptotic given by Eq.ll. 

Another special value of the mechanical energy is 
E(r) = 0, corresponding to parabolic motion. In this case, 



(11) y par = tf A i?[i + 2 (^) 3 ] 1/2 



R 



(14) 



The two special trajectories with E{r) = E m i n and E{r) = 
(see Fig.l) converge to the asymptotic of Eq.ll at large 
distances - in agreement with the results of Sec. 3. 

It is quite interesting that the most accurate data on the 
flow (from the TRGB and Cepheids; dots in Fig.l) show a 
low dispersion around the line V = HyR, and the cor- 
responding trajectories occupy the area just between the 
curves Vm n (i?) and V^ ar (i?). The rather wide spread of the 
other data in Fig.l is mostly due to observational errors 
(Karachentsev & Nasonova 2010 - especially their Fig. 7b). 
Anyway, the restriction E m i n < E < puts an upper limit 
on the value \E\ in Eq.5, and because of this, the asymptotic 
of Eq.ll may be reached more easily in the flow. 

It is also suggestive that the Local Group outflow 
has the same feature: the galaxies prefer the area be- 
tween Vmin and V par in the Hubble diagram of Fig. 2 (c.f. 
Teerikorpi et al. 2008). This is another similarity aspect 
between the Virgo system and the Local system (Sec. 2). 
Mathematically, the similarity reflects the fact that Eqs.13 
and 14 have exactly the same view for both systems in the 
dimensionless form y(x) = with y = V/(H a Rzg), % = 
R/Rzg', the scale factor is again about 10. 

This similarity suggests that the condition i? m i n < E < 
has more general significance; it may reflect the kinemat- 
ical state of the flow at the early formation epoch. For ex- 
ample, V > Vmin follows, if the galaxies of the flow escaped 
initially from the cluster (Chernin et al. 2004, 2007c). 



4.3. Comparison with cosmology 

The zero-gravity radius Rzg is a local spatial counterpart of 
the redshift z A ~ 0.7 at which the gravity and the antigrav- 
ity balance each other in the Universe as a whole. Globally, 
the balance takes place only at one moment t(z A ) in the 
entire co-moving space, while the local gravity-antigravity 
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balance exists during all life-time of the system since its 
formation, but only on the sphere of the radius Rzg- 

The Virgocentric flow should not be identified with the 
global cosmological expansion. The Virgo system lies deeply 
inside the cosmic cell of uniformity and does not "know" 
about the universe of the horizon scales. In the early epoch 
of weak protogalactic perturbations, the dark matter and 
baryons of the system participated in the cosmological ex- 
pansion. But later, at the epoch of nonlinear perturba- 
tions, the matter separated from the expansion and under- 
went violent evolution. Contrary to that, the unperturbed 
flow at horizon-scale distances was not similarly affected 
and up to now the initial isotropy of its expansive motion 
and the uniformity of its matter distribution are conserved. 
Therefore, it is mysterious that the global Hubble constant 
72 km/s/Mpc is so close to the local expansion rate 65 
km/s/Mpc in the Virgocentric flow (two beams in Fig.l). 
This is one aspect of the Hubble-Sandage paradox (Sec.l). 

Trying to find an explanation, we note that the 
Friedmann equation for the scale factor a(t) has much the 
same form as the local energy conservation equation (5): 

\a{tf = ^L )+ ^G PA a{tf+B. (15) 

The constant C = ^-pM(t)a(t) 3 has the dimension of mass, 
and pm (t) is the density of matter (dark and baryonic) ; the 
energy constant B — in the standard flat cosmology. The 
cosmological Hubble factor comes from this equation: 

H{t) = m. = [^ (pmW+pa)] i/2 = HA[l+ ^^ m 

When time goes to infinity, the matter density drops to 
zero, the dark energy becomes dominant, and the cosmo- 
logical Hubble factor becomes time-independent: 

ff(t)^F A = [^-p A ] 1/2 = 63-65 km/s/Mpc. (17) 

At the present epoch, tu = 13.7 Gyr, the Hubble rate 

H = H A [l + ^ L ] 1/2 -72 km/s/Mpc, (18) 

"A 

which is equal to the empirical (WMAP) value. 

As we see, the Virgocentric flow and the global cosmo- 
logical expansion have a common asymptotic in the limit 
of dark energy domination: for the local flow, this is the 
asymptotic in space (large distances from the cluster cen- 
ter), and for the global flow, this is the asymptotic in time 
(a future of the Universe). If the asymptotic is approached 
now both globally and locally, we would have H ~ H\ on 
the largest scales and Hy ~ Ha in the Virgocentric flow. 

In their present states, both flows are near the asymp- 
totic, though do not quite reach it. That is why H, Hy, and 
Ha have similar values. The flows are not related to each 
other directly; however, they "feel" the same uniform dark 
energy background, and in this way both are controlled by 
the antigravity force which is equally effective near the cos- 
mic horizon and at the periphery of the Virgocentric flow. 
As a result, the local flow at its large distances looks much 
like a "part" of the global horizon-distance flow. 

The near coincidence of the values H, Hy, and Ha ex- 
plained by our model clarifies the Hubble-Sandage paradox 
(Sec.l) for the Virgocentric flow: asymptotically, the flow 
reaches not only the linear Hubble law (Sec. 3), but also ac- 
quires the expansion rate Hy which approaches the global 
rate when the measuring accuracy is improved. 



4.4. Probing local dark energy 

The global density of dark energy can be derived, if the 
redshift za is found at which the gravity of matter is ex- 
actly balanced by the dark energy antigravity, and also the 
present matter density is known. The same logic works in 
local studies: just find the zero-gravity radius Rzg from the 
outflow observations and the mass of the local system. 

We may now restrict the value of Rzg using the diagram 
of Fig.l. Since the zero-gravity surface lies outside the clus- 
ter volume, it should be that Rzg > Ro — 6 — 7 Mpc. On 
the other hand, the linear velocity-distance relation, with 
the Hubble ratio close to H, is clearly seen from a distance 
of about 15 Mpc. This suggests that Rzg < 15 Mpc. This 
argument that the zero-gravity surface is located below the 
point where the local flow reaches the global expansion 
rate, gained support from the calculations by Teerikorpi 
& Chernin (2010). As for the cluster mass, we adopt the 
range M = (0.6 - 1.2) x 10 15 M Q (Sec.2). Then Eq.6 in the 
form p x — M/(^-R^ G ) directly leads to robust upper and 
lower limits to the unknown local dark energy density p x : 

0.2<p I <5x 10~ 29 g/cm 3 . (19) 

In fact, the upper limit obtained from the cluster size (Rzg) 
is conservative. We note that the theoretical minimum ve- 
locity curve in Fig.l would be shifted up to the velocity- 
distance relation defined by the TRGB and Cepheid dis- 
tances, if one takes p x = 1.4 x 10~ 29 g/cm 3 . This would be 
a more realistic upper limit, and near the limit similarly ob- 
tained in Teerikorpi et al. (2008) using the Local Group and 
the M81 group. Thus the local dark energy density is near 
the global value or may be the same; pa = 0.73 x 10~ 29 
g/cm 3 lies comfortably in the range of Eq.19. Similar es- 
timates for the dark energy come from studies of the 
Local Group, M81 group and Cen A group (Chernin et 
al. 2007a,b, 2009). 

We may now argue the other way round: let us assume 
that the density of dark energy is known from cosmological 
observations and the zero-gravity radius Rzg is between 6 
and 15 Mpc. Then we may estimate independently the mass 
of the Virgo Cluster: M = (0.3 - 4) x 10 15 M Q . The value 
of the cluster mass adopted earlier is within this interval. 

4.5. The origin of the local flows? 

The model of Sec. 3 describes the present structure of the 
Virgocentric flow and also predicts its future evolution. 
However, it says almost nothing about the initial state of 
the flow. Our discussion suggests that the Virgo cluster 
and the outflow form a system with a common origin and 
evolutionary history. Its formation was most probably due 
to complex linear and nonlinear processes (collisions and 
merging of galaxies and protogalactic units, violent relax- 
ation of the system in its self-gravity field, etc.). Similar pro- 
cesses are expected in the early history of the Local Group 
and other systems with the two-component structure. Their 
basic features might be recognized from big N-body simu- 
lations, like the Millennium Simulation (Li & White 2008). 
We leave it for a later study to clarify the way in which the 
galaxies of local flows gain their initial expansion velocities. 
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5. Conclusions 

The recently published systematic and most accurate data 
on the velocities and distances of galaxies in the Virgo clus- 
ter and its environment (Karachentsev & Nasonova 2010) 
shed light on the physics behind the observed properties of 
the cluster and the flow. We find that: 

1. The cluster and the flow can be treated together as a 
physical system embedded in the uniform dark energy back- 
ground; this is adopted in our analytical nonlinear model. 

2. The nonlinear interplay between the gravity produced 
by the cluster mass (dark matter + baryons) and the anti- 
gravity of the dark energy is the major dynamical factor 
determining the kinematic structure of the system and con- 
trolling its evolution; our model enables one to describe this 
in a simple (exact or approximate) quantitative way. 

3. The key physical parameter of the system is its zero- 
gravity radius, Rzg = 9 — 11 Mpc, below which the grav- 
ity dominates the quasi-stationary bound cluster, while the 
antigravity dominates the expanding Virgocentric flow. 

4. The dark energy antigravity brings order and regu- 
larity to the Virgocentric flow and the flow acquires the lin- 
ear velocity-distance relation at the Virgocentric distances 
R > 15 Mpc. Thus, the Hubble-Sandage paradox is under- 
stood on the cluster scale. 

5. In the velocity-distance diagram, the zero-gravity ra- 
dius is near to or somewhat larger than the zero-velocity 
radius Rq = 5.0 — 7.5 Mpc, as determined by Karachentsev 
& Nasonova (2010). 

6. With the value of Rzg <1 Ro found from the di- 
agram, one may estimate the mass of the Virgo cluster, 
using the known global density of dark energy: M = 
(0.3-4) x 1O 15 M . 

7. With the same value of Rzg, one may estimate the 
local density of the dark energy in the Virgo system, if the 
mass of the cluster is considered known (from the virial or 
the zero-velocity method): p\ — (0.2 — 5) x 10 _29 g/cm 3 . 

8. The two-component Virgo system is similar to the 
Local system formed by the Local Group together with the 
local outflow around it. In the velocity-distance diagram, 
the Virgo system reproduces the Local system with the 
scale factor of about 10. A dynamical similarity is seen in 
the fact that the zero-gravity radius for the Virgo system 
is about 10 times larger than for the Local system. 

9. The similarity suggests that a universal two- 
component grand design may exist with a quasi-stationary 
bound central component and an expanding outflow around 
it on the spatial scales of 1-30 Mpc. The phase and dynam- 
ical structure of the real two-component systems reflects 
the nonlinear gravity-antigravity interplay with dark en- 
ergy domination in the flow component. 
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